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£f~^ , The Muon g — 2 collaboration has measured the anomalous magnetic g value, a = (g — 2)/2, 

of the positive muon with an unprecedented uncertainty of 0.7 parts per million. The result 
• a fl+ (expt) = 11659 204(7)(5) x 10~ 10 , based on data collected in the year 2000 at Brookhaven 

t-H \ National Laboratory, is in good agreement with the preceding data on and a^- . The measure- 

ment tests standard model theory, which at the level of the current experimental uncertainty involves 

C^) ' quantum electrodynamics, quantum chromodynamics, and electroweak interaction in a significant 

C ~ ) ' way. 

I. INTRODUCTION 

V ■ ... .... 

^ c"| The magnetic moment fl of a particle with charge e, mass m, and spin s is given by 

' e 

• £h ! P= o (1) 

2mc 

in which g is the gyromagnetic ratio. For point particles with spin-i, Dirac theory predicts g = 2. 

Precision measurements of the g values of leptons and baryons have historically played an important role in the 
development of particle theory. The g value of the proton, for example, is found to differ sizeably from 2, which provides 
evidence for a rich internal proton (spin) structure. The lepton g values deviate from 2 only by about one part in a 
thousand, consistent with current evidence that leptons are point particles with spin-i. The anomalous magnetic g 
value of the electron, a e — (g e — 2)/2, is among the most accurately measured quantities in physics, and is presently 
known with an uncertainty of about four parts per billion (ppb) . The value of a e is described in terms of standard 
model (SM) field interactions, in which it has a leading order contribution of a/(27r) from the so-called Schwinger 
term. Nearly all of the measured value is contributed by QED processes involving virtual electrons, positrons, and 
photons. Particles more massive than the electron contribute only at the level of the present experimental uncertainty. 

The anomalous magnetic g value of the muon, a^, is more sensitive than a e to processes involving particles more 
massive than the electron, typically by a factor (m^/rrie) 2 ~ 4 • 10 4 . A series of three experiments at CERN 
measured to a final uncertainty of about 7 parts per million (ppm) , which is predominantly of statistical origin. The 
CERN generation of experiments thus tested electron- muon universality and established the existence of a ~ 59 ppm 
hadronic contribution to a^. Electroweak processes are expected to contribute to at the level of 1.3 ppm, as are 
many speculative extensions of the SM. 

The present muon g — 2 experiment at Brookhaven National Laboratory (BNL) has determined of the positive 
muon with an uncertainty of 0.7 ppm from data collected in the year 2000, and aims for a similar uncertainty on a M - 
of the negative muon from measurements in 2001. The continuation of the experiment by a single running period, if 
funded, should bring the design goal precision of 0.4 ppm within reach. 



FIG. 1: Top view of the g — 2 apparatus. The beam of longitudinally polarized muons enters the superferric storage ring magnet 
through a superconducting inflector magnet located at 9 o'clock and circulates clockwise after being placed onto stored orbit 
with three pulsed kickers modules in the 12 o'clock region. Twenty-four lead scintillating-fiber calorimeters on the inner, open 
side of the C-shaped ring magnet are used to measure muon decay positrons. The central platform supports the power supplies 
for the four electrostatic quadrupoles and the kicker modules. 



II. EXPERIMENT 

The concept of the experiment at BNL is the same as that of the last of the CERN experiments [H-d-tH and involves 
the study of the orbital and spin motions of polarized muons in a magnetic storage ring. 

The present experiment (Fig. ^) is situated at the Alternating Gradient Synchrotron (AGS), which in the year 
2000 delivered up to 60 x 10 12 protons in twelve 50 ns (FWHM) bunches over its 3 s cycle. The 24GeV protons from 
the AGS were directed onto a rotating, water-cooled nickel target. Pions with energies of 3.1 GeV emitted from the 
target were captured into a 72 m straight section of focusing-defocusing magnetic quadrupoles, which transported 
the parent beam and naturally polarized muons from forward pion decays. At the end of the straight section, the 
beam was momentum-selected and injected into the 14.2 m diameter storage ring magnet [f| through a field-free 
inflector region in the magnet yoke. A pulsed magnetic kicker Q located at approximately one quarter turn from 
the inflector region produced a lOmrad deflection which placed the muons onto stored orbits. Pulsed electrostatic 
quadrupoles [8j provided vertical focusing. The magnetic dipole field of about 1.45 T was measured with an NMR 
system |!j relative to the free proton NMR frequency uj p over most of the 9 cm diameter circular storage aperture. 
Twenty-four electromagnetic calorimeters [l(j read out by 400 MHz custom waveform digitizers (WFD) were used on 
the open, inner side of the C-shaped ring magnet to measure muon decay positrons. The decay violates parity, which 
leads to a relation between the muon spin direction and the positron energy spectrum in the laboratory frame. For 
positrons above an energy threshold E, the muon-decay time-spectrum 

N(t) = N (E) exp (^) [1 + A(E) sin {uj a t + <j>(E))\ , (2) 

in which Nq is a normalization, 7T ~ 64 /zs is the dilated muon lifetime, A an asymmetry factor, a phase, and uj a 
the angular difference frequency of muon spin precession and momentum rotation. In our measurements, the NMR 
and WFD clocks were phase-locked to the same LORAN-C 01 frequency signal. 

The muon anomalous magnetic g value is evaluated from the ratio of the measured frequencies, R = Lu a /uj p , 
according to: 

a, = ^ (3) 

in which A = fj,^/ \x v is the ratio the muon and proton magnetic moments. The value with smallest stated uncertainty, 
A = = 3.183 345 39(10) [T^ | , results from measurements of the microwave spectrum of ground state muonium [l3| 

and theory Q [l5| . 
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TABLE I: Systematic uncertainties for the cu p analysis. The uncertainty "Others" groups uncertainties caused by higher 
multipoles, the trolley frequency, temperature, and voltage response, eddy currents from the kickers, and time-varying stray 
fields. 



Source of errors 


Size [ppm] 


Absolute calibration of standard probe 


0.05 


Calibration of trolley probe 


0.15 


Trolley measurements of Bo 


0.10 


Interpolation with fixed probes 


0.10 


Uncertainty from muon distribution 


0.03 


Others 


0.10 


Total systematic error on uj p 


0.24 



Important improvements made since our preceding measurement [16j include: the operation of the AGS with 12 
beam bunches, which contributed to a 4-fold increase in the data collected; a new superconducting inflector magnet, 
which greatly improved the homogeneity of the magnetic field in the muon storage region; a sweeper magnet in 
the beamline, which reduced AGS background; additional muon loss detectors, which improved the study of time 
dependence; and further refined analyses, in particular of coherent betatron oscillations. 

III. DATA ANALYSIS 

The analysis of a M follows, naturally, the separation of the measurement in the frequencies uj p and oj a . Both 
frequencies were analyzed independently by several groups within the collaboration. The magnetic field frequencies 
measured during the running period were weighted by the distribution of analyzed muons, both in time and over 
the storage region. The frequency fitted from the positron time spectra was corrected by +0.76(3) ppm for the net 
contribution to the muon spin precession and momentum rotation caused by vertical beam oscillations and, for muons 
with 7 ^ 29.3, by horizontal electric fields ^tJ- The values of R — uj a /ui p and a M were evaluated only after each of 
the frequency analyses had been finalized; at no earlier stage were the absolute values of both frequencies, ui p and u a , 
known to any of the collaborators. 

A. The frequency ui p 

The analysis of the magnetic field data starts with the calibration of the 17 NMR probes in the field trolley using 
dedicated measurements taken during and at the end of the data collection period. In these calibration measurements, 
the field in the storage region was tuned to very good homogeneity at two specific calibration locations. The field 
was then measured with the NMR probes mounted in the trolley shell, as well as with a single probe plunged into the 
storage vacuum and positioned to measure the field values in the corresponding locations. Drifts of the field during 
the calibration measurements were determined by remeasuring the field with the trolley after the measurements with 
the plunging probe were completed, and in addition by interpolation of the readings from nearby NMR probes in the 
outer top and bottom walls of the vacuum chamber. The difference of the trolley and plunging probe readings forms 
a calibration of the trolley probes with respect to the plunging probe, and hence with respect to each other. The 
plunging probe, as well as a subset of the trolley probes, were calibrated with respect to a standard probe at the 
end of the running period in a similar sequence of measurements in the storage region, which was opened to air for 
that purpose. The leading uncertainties in the calibration procedure result from the residual inhomogeneity of the 
field at the calibration locations, and from position uncertainties in the active volumes of the NMR probes. These 
uncertainties were evaluated from measurements in which the trolley shell was purposely displaced and known field 
gradients were applied using the so-called surface and dipole correction coils of the ring magnet. The size of these 
uncertainties is estimated to be 0.15 ppm, as listed in Table |U The uncertainty in the absolute calibration of the 
standard probe amounts to 0.05 ppm [Ti|. The dependencies of the trolley NMR readings on the supply voltage and 
on other parameters were measured to be small in the range of operation. An uncertainty of 0.10 ppm ("Others" in 
Table Pi is assigned, which includes also the measured effects from the transient kicker field caused by eddy currents 
and from AGS stray fields. 

The magnetic field inside the storage region was measured 22 times with the field trolley during the data collection 
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Multipoles [ppm] 
normal skew 
Quad 0.24 0.29 
Sext -0.53 -1.06 
Octu -0.10 -0.15 
Decu 0.82 0.54 



azimuth [degree] radial distance [cm] 



FIG. 2: The NMR frequency measured with the center trolley probe relative to a 61.74 MHz reference versus the azimuthal 
position in the storage ring (left), and (right) a 2-dimensional multipole expansion of theazimuthal average of the field measured 
with 15 trolley probes with respect to the central field value of 1.451 275 T. The multipole amplitudes are given at the storage 
ring aperture, which has a 4.5 cm radius as indicated by the circle. 



from January to March 2000. Fig. |3l shows the field value measured in the storage ring with the center trolley 
probe versus the azimuthal angle. The field is seen to be uniform to within about ±50 ppm of its average value 
over the full azimuthal range, including the region near 350° where the inflector magnet is located. Non-linearities 
in the determination of the trolley position during the measurements — from the measured cable lengths and from 
perturbations on the readings from fixed probes as the trolley passes — are estimated to affect the azimuthal average 
of the field at the level of 0.10 ppm. Fig. [2b shows a 2-dimensional multipole expansion of the azimuthally averaged 
readings from the trolley probes, 

oc oo 

B y = J2 C nr n cos(ncl>)- J2 D nr n sm(n<f>), (4) 
n=0 n=0 

oo oo 

B x = ^<7 n r"sin(n^)+^D„r"cos(n0), (5) 

n=0 ra=0 

where the coefficients C n and D n are the normal and skew multipoles, and r and <f> denote the polar coordinates in the 
storage region. The multipole expansion was truncated in the analysis after the decupoles. Measurements with probes 
extending to larger radii show that the neglect of higher multipoles is at most 0.03 ppm in terms of the average field 
encountered by the stored muons, in agreement with magnet design calculations. The field averaged over azimuth is 
seen to be uniform to within 1.5 ppm of its value. 

The measurements with the trolley relate the readings of 370 NMR fixed probes in the outer top and bottom walls 
of the storage vacuum chamber to the field values in the beam region. The fixed NMR probes are used to interpolate 
the field when the field trolley is 'parked' in the storage vacuum just outside the beam region, and muons circulate in 
the storage ring. Since the relationship between the field value in the storage region and the fixed probe readings may 
change during the course of the data collection period, the field mappings with the trolley were repeated typically two 
to three times per week, and whenever ramping of the magnet or a change in settings required such. The uncertainty 
associated with the interpolation of the magnetic field between trolley measurements is estimated from the spread of 
the difference between the dipole moments evaluated from the fixed probe measurements and from the trolley probe 
measurements in periods of constant magnet settings and powering. It is found to be 0.10 ppm. 

Since the field is highly uniform, the field integral encountered by the (analyzed) muons is rather insensitive to the 
exact location of the beam. As in earlier works |!6lll9| . the radial equilibrium beam position was determined from the 
debunching of the beam following injection and the vertical position from the distribution of counts in scintillation 
counters mounted on the front faces of the positron calorimeters. The position uncertainty amounts to 1 - 2 mm, 
which contributes 0.03 ppm uncertainty to the field integral. 

The result for field frequency oj p weighted by the muon distribution is found to be, 

lu p /(2tt) = 61 791 595(15) Hz (0.2 ppm), (6) 

where the uncertainty has a leading contribution from the calibration of the trolley probes and is thus predominantly 
systematic. A second analysis of the field has been performed using additional calibration data, a different selection 
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FIG. 3: The time spectrum for 4 • 10 positrons with energies greater than 2 GeV collected from January to March 2000, 
after corrections for pile-up and for the bunched time structure of the injected beam (left) were made, and (right) the Fourier 
transform of the time spectrum, in which muon decay and spin precession (cf. Eq. 2) has been suppressed to emphasize other 
effects. 

of fixed NMR probes, and a different method to relate the trolley and fixed probe readings. The results from these 
analyses are found to agree to within a fraction of the total uncertainty on w p . 



B. The frequency ui a 



The event sample available for analysis from data collection in the year 2000 amounts to about 4 • 10 9 positrons 
reconstructed with energies greater than 2 GeV and times between 50 fis and 600 fis following the injection of a beam 
bunch. Fig. shows the time spectrum corrected for the bunched time structure of the beam and for overlapping 
calorimeter pulses, so called pile-up [16J. 

The leading characteristics of the time spectrum are those of muon decay and spin precession (cf. Eq.[5J). Additional 
effects exist, as seen from the Fourier spectrum in Fig. [3Jd, and require careful consideration in the analysis. These 
effects include detector gain and time instability, muon losses, and oscillations of the beam as a whole, so-called 
coherent betatron oscillations (CBO). 

Numerically most relevant to the determination of Lo a are CBO in the horizontal plane. CBO are caused by injecting 
the beam through the relatively narrow 18(w) x 57(h) mm 2 aperture of the 1.7 m long inflector channel into the 90 mm 
diameter aperture of the storage region, and have been observed directly with fiber harp monitors plunged into the 
beam region for this purpose. The CBO frequency is determined by the focusing index of the storage ring, and is 
numerically close to twice the frequency uj a for the quadrupole settings employed in most measurements so far. Since 
the calorimeter acceptances vary with the radial muon decay position in the storage ring and with the momentum of 
the decay positron, the time and energy spectra of the observed positrions are modulated with the CBO frequency. 
These modulations affect the normalization No, the asymmetry A, and the phase 4> in Eq. |2Jat the level of 1%, 0.1%, 
and 1 mrad at beam injection. When not accounted for in the function fitted to the data, the modulations of the 
asymmetry and phase with a frequency w c bo,h - 2xw a may manifest themselves as artificial shifts of up to 4 ppm 
in the frequency values lu u determined from individual calorimeter spectra. The circular symmetry of the experiment 
design results in a strong cancellation of such shifts in the joined calorimeter spectrum. 

Several approaches have been pursued in the analysis of u> a . In one approach, the time spectra from individual 
positron calorimeters was fitted in narrow energy intervals using a fit function as in Eq. [21 extended by the number, 
asymmetry, and phase modulations. Other approaches made use of the cancellation in the joined calorimeter spectra 
and either fitted for the residual of the leading effects, or accounted for their neglect in a contribution to the systematic 
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TABLE II: Systematic uncertainties for the ui a analysis. The uncertainty "Others" groups uncertainties caused by AGS 
background, timing shifts, vertical oscillations and radial electric fields, and beam debunching/randomization. 



Source of errors 


Size [ppm] 


Coherent betatron oscillations 


0.21 


Pileup 


0.13 


Gain changes 


0.13 


Lost muons 


0.10 


Binning and fitting procedure 


0.06 


Others 


0.06 


Total systematic error on u a 


0.31 



uncertainty. The results are found to agree, on u a to within the expected 0.5 ppm statistical variation resulting from 
the slightly different selection and treatment of the data in the respective analyses. The combined result is found to 
be, 

w /(2tt) = 229074 11(14)(7) Hz (0.7 ppm), (7) 

in which the first uncertainty is statistical and the second systematic. The combined systematic uncertainty is broken 
down by source in Table ITT1 



IV. RESULTS AND DISCUSSION 



The value of a M was evaluated after the analyses of to p and oj a had been finalized, 

= 11659 204(7)(5) x 10~ 10 (0.7ppm), (8) 

where the first uncertainty is statistical and the second systematic. This new result is in good agreement with the 
previous measurements 0, [lg, Ha, 120] and drives the present world average, 

a M (exp) = 11659 203(8) x 10" 10 (0.7ppm), (9) 

in which the uncertainty accounts for known correlations between the systematic uncertainties in the measurements. 
Fig. ^] shows our recent measurements of , together with two SM evaluations discussed below. 

In the SM, the value of a M receives contributions from QED, hadronic, and electroweak processes, a M (SM) = 
a /J (QED) + (had) + a M (weak). The QED and weak contributions can, unlike the hadronic contribution, be evaluated 
perturbatively, a Al (QED) = 11658 470.57(29) x 10~ 10 ^ and a M (weak) = 15.1(4) x 10" 10 The hadronic 

contribution is, in lowest order, related by dispersion theory to the hadron production cross sections measured in 
e + e~ collisions and, under additional assumptions, to hadronic r-decay. Clearly, the hadronic contribution has a long 
history of values as new data appeared and analyses were refined. 

Shortly before the SPIN-2002 conference, Davier and co-workers released a new and detailed evaluation |2l|, which 
now incorporates the high precision e + e~ data [24j in the region of the p resonance from CMD-2 at Novosibirsk, more 
accurate e + e~ measurements |25l l26j in the 2-5 GeV energy region from BES in Beijing, preliminary results from 
the final ALEPH analysis [^3 of hadronic r-decay at LEP1, as well as additional CLEO data [2^,|2!|. The authors 
note discrepancies between the e + e~ and r data at the present levels of precision, and obtain separate predictions 
for the contribution to a p (SM) from lowest order hadronic vacuum polarization, a p (had, 1) = 685(7) x 10~ 10 from 
e + e _ data and a,,, (had, 1 ) = 702(6) x 10~ 10 from r data. Higher order contributions include higher order hadronic 
vacuum polarization [30, El and hadronic light-by- light scattering [3^, |U 0, HE US • 

Open questions concern the SM value of a M , in particular the hadronic contribution, and the experimental value 
of a„- at sub-ppm precision. The former should benefit from further theoretical scrutiny, e.g. Refs. |37l l38t l39l [40j. 
from radiative-return measurements at e + e~ factories, e.g. Refs. ^3,E2, an( A possibly from latice calculation |43| . 
We are currently analyzing a sample of about 3 x 10 9 decay electrons from a„- data collected in the year 2001. Stay 
tuned! 
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FIG. 4: Recent measuremens of a M and standard model evaluations using the evaluations in Ref. 
contribution from hadronic vacuum polarization. 
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